Violation of energy-per-hadron scaling in a resonance matter 



L. Bravina, 1 '^ Amand Faessler, 1 C. Fuchs, 1 Zhong-Dao Lu, 2,3 and E.E. Zabrodin 1 '^ 

1 Institut fur Theoretische Physik, Universitat Tubingen, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 

2 China Institute of Atomic Energy, P. O. Box 275(18), Beijing 1024-13, China 

3 China Center of Advance Science and Technology (CCAST), Beijing 100080, China 

(Dated: February 8, 2008) 

Yields of hadrons, their average masses and energies per hadron at the stage of chemical freeze-out 
in (ultra)relativistic heavy-ion collisions are analyzed within the statistical model. The violation of 
the scaling {E)/(N) = 1 GeV observed in Au+Au collisions at y/s = 130 A GeV is linked to the 
formation of resonance-rich matter with a considerable fraction of baryons and antibaryons. The 
rise of the energy-per-hadron ratio in baryon-dominated matter is discussed. A violation of the 
scaling condition is predicted for a very central zone of heavy-ion collisions at energies around 40^4 
GeV. 
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The properties of nuclear matter under extreme con- 
ditions have been the subject of intensive experimental 
and theoretical studies during the last few decades. Up 
to now experiments with heavy-ion collisions remain the 
only means to explore the properties of hot and dense nu- 
clear matter in the laboratory. Two energy ranges con- 
nected with predicted phase transitions in the nuclear 
medium have been studied especially vigorously, - the 
intermediate energy range, where the nuclear liquid-gas 
phase transition might occur, and the range of relativistic 
and ultra-relativistic energies, where the phase transition 
to a deconfined quark-gluon plasma (QGP) and a restora- 
tion of chiral symmetry should take place (see jl], ^ and 
references therein). The search for the QGP formation 
is one of the top-priority goals of the heavy-ion collider 
program at Relativistic Heavy Ion Collider (RHIC) in 
Brookhavcn, which is operating since 1999, and at the 
forthcoming Large Hadron Collider (LHC) at CERN. 

In order to reveal possible fingerprints of the QGP for- 
mation, various microscopic (transport, string, cascade) 
and macroscopic (thermal and hydrodynamic) models 
are widely used for the analysis of measured particle 
abundances and energy spectra. In macroscopic scenar- 
ios it is assumed that a rapidly expanding and cooling 
thermalized system experiences at the late stage of its 
evolution the so-called chemical freeze-out, where all in- 
elastic processes have to cease, accompanied by the ther- 
mal freeze-out, which occurs when the mean free path of 
particles exceeds the linear sizes of the system. There- 
fore, the conditions of the system at the chemical freeze- 
out stage can be obtained from hadron abundances and 
ratios, which are not affected by the collective flow. The 
analysis of experimental data taken in a broad energy 
range from SIS to SPS suggests || that there exists a 
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scaling law concerning the energy per hadron ratio at the 
chemical freeze-out, (E)/(N) ps 1 GeV, while in Au+Au 
collisions at y/s — 130 ^4GeV (RHIC) this ratio increases 
to (E) I (N) ps 1.1 GeV. Simple estimates show that in 
the latter case the average hadron mass should increase 
as well, because the chemical freeze-out temperature at 
RHIC does not exceed that at SPS by more than 10- 
20 MeV 0. In the present paper we argue that the 
rise of the average hadron energy and mass at RHIC is 
caused by the transition to dense meson-resonance rich 
matter with a considerable fraction of baryons and an- 
tibaryons. Another substance with peculiar character- 
istics is baryon-dominated resonance matter, which can 
be formed in heavy-ion collisions at bombarding energies 
between 10 ^4GeV and 40 j4GeV, accessible for the accel- 
erator planned at GSI. Here the energy per hadron at 
chemical freeze-out increases as well. 

For our study we use a conventional statistical model 
(SM) of an ideal hadron gas g, |, |, |, @, § § |fij 
which enables one to determine all macroscopic char- 
acteristics of a system at given temperature T, baryon 
chemical potential fis, and strangeness chemical poten- 
tial us via the set of distribution functions (in units of 
c = ks = h = 1) 
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Here to;, Ei — \Jp 2 + mf, p, and /i; are the mass, energy, 
momentum, and the chemical potential of hadron species 
i, respectively. Denoting baryon and strange charges of 
the ith particles as Bi and Si one can write fii = B^b + 
Sifis- The electrochemical potential considered in [lC ] 
and the isospin chemical potential considered in 
are neglected. The sign "+(—)" in Eq. ([!]) corresponds 
to fermions (bosons). The particle number density n^, 
the energy density £; and the partial pressure Pi read 
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p * = WrJ wt^ mmi)dp > (4) 

with gi being the spin-isospin degeneracy factor of hadron 
i. Instead of evaluating of the integrals in Eqs. 
by replacing of Fermi-Dirac or Bose-Einstein distribution 
functions (El) to Maxwell-Boltzmann ones 
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we employ the series expansion of (Q) in a form 
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which is inserted to Eqs. (|])-([|). After some straightfor- 
ward calculations one gets 
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where iiTi and K% are modified Hankel function of first 
and second order. For n = 1 the results with the 
Maxwell-Boltzmann distribution (||) are regained |p^| . 

The equation of state (EOS) of an ideal hadron gas in 
the form (E(T, hb)) /(N(T, ^b)) and the average hadron 

N 

mass (M(T, hb)) — l/N^^rm are shown in Fig. |l|. The 

calculations are performed within the standard SM at 
zero net strangeness in the system. The increase of both 
quantities can be caused either by the rise of the temper- 
ature of the system or by an increasing baryon chemical 
potential, which is directly linked to the baryon density. 
Both macroscopic [Q] and microscopic |l3| models indi- 
cate that the case with low \xb and relatively high T is 
relevant for Au+Au collisions at RHIC. To study the con- 
ditions of chemical freeze-out we also apply a generaliza- 
tion of the SM to two thermal sources (TSM) proposed in 
flif , which accounts for possible inhomogeneities of the 
net baryon charge, observed experimentally e.g. at SPS 
energies (158 AGeV/c) fll5fl , and net strangeness distri- 
bution inside the reaction volume. The idea behind the 
TSM is quite simple. It is well known (see, e.g., ||]) 
that the hadron ratios in the system consisting of several 
fireballs are not affected by longitudinal or transverse 
collective motion of sources, provided all fireballs have 




FIG. 1: (a) Average energy per hadron and (b) average 
hadron mass as functions of the temperature T and the baryo- 
chemical potential \ib of the system. 



the same baryon chemical potential and the same tem- 
perature. In this case the particle ratios are identical to 
those obtained for a single static source. If, for instance, 
the baryon charge or strangeness is not homogeneously 
distributed within the total volume, the problem cannot 
be reduced to a single source scenario. In the TSM the 
whole system is separated in two parts, which at SPS 
and lower energies can most naturally be interpreted as 
an inner source (or core) and an outer source (or halo), 
each being in local thermal and chemical equilibrium, i.e. 
their (local) macroscopic characteristics are determined 
by Eqs. (||) - (||). However, their temperatures, baryon 
and strangeness chemical potentials are allowed (but not 
postulated) to be different. No additional constraints 
are assumed in the model except of the total strangeness 
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conservation Ng ot — 

A fit to experimental data at SPS and RHIC ener- 
gies has been performed in [|l4|. Here it turns out that 
at RHIC the results of the fit to the SM and the TSM 
are identical: the TSM simply splits the volume of the 
system in two equal parts with similar characteristics, 
which means that the particles are really emitted from 
one thermalized and homogeneous source. This gives us 
a temperature of T = 176 MeV and a baryo-chemical po- 
tential of fiB = 39.8 MeV |I6| |. The energy per hadron 
and the average hadron mass are (E) / (N) = 1.12 and 
(AI) = 0.77 GeV, respectively. As seen in Fig. 0, this 
is a significant rise compared to the values of (E) / (N) 
and (M) obtained from the fit to experimental data on 
heavy- ion collisions at lower energies. 

In contrast to RHIC, at lower energies results of the 
SM and the TSM do not coincide anymore, although both 
models seem to provide reasonable agreement with the 
data, as shown in Table ||. Here it is worth noting that 
the TSM analysis has been performed for RHIC data 
taken in a very narrow midrapidity window. However, 
results of the TSM fitted separately to merely 47r-data 
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FIG. 2: (a) Average energy per hadron and (b) average 
hadron mass as functions of the center-of-mass energy of col- 
liding nuclei. Full circles correspond to central source of the 
TSM, open circles denote the SM results. Lines are drawn to 
guide the eye. 



and to midrapidity data, obtained in Pb+Pb collisions 
at SPS energies, indicate jl4j that the temperature and 
the baryon chemical potential of the central source are 
varying insignificantly. Therefore, we do not expect con- 
siderable changes of the core conditions when the 47r-data 
at RHIC will be available. Macroscopic characteristics 
of the system obtained from the fit to the models are 
listed in Table 0. For A+A collisions at E lab = 11.6 
AGeV, 40AGeV, and 80AGeV the energy per hadron 
and the average hadron mass in the central source (using 
the TSM) are about 20% larger than those obtained in 
the SM because of the following reasons: At such ener- 
gies the TSM favours the formation of hot and dense core 
surrounded by a cooler halo. Moreover, the model indi- 
cates that the net baryon charge and the net strangeness 
are non-uniformly distributed within the system volume. 
Here three points should be mentioned: Using the TSM, 
all antibaryons are contained in the core, which is in ac- 
cord with experimental data on, e.g., p, d, and A pro- 
duction |l^, 2^, ^tJ. The net strangeness in the core 
is small and negative. Such a scenario is confirmed by 
microscopic ultra-relativistic quantum molecular dynam- 
ics (UrQMD) m model predictions (29). The negative 
strangeness in the central zone of a heavy-ion collision at 
energies below RHIC can be explained by different inter- 
action cross sections of kaons and antikaons with baryons. 
In contrast, at RHIC energies and above the medium is 
meson dominated and the antibaryon yield at midrapid- 
ity is close to that of baryons. Hence, differences in the 



TABLE I: Hadron yields and ratios for central heavy-ion 
collisions at 11.64 GeV, 404 GeV, and 804 GeV, respectively, 
and results of the fit to the two-source and the single-source 
statistical models of an ideal hadron gas. 



Data 



TSM 



SM Ref. 



Au4 
N B 

p/lT' 
7T+ 

K+ 
K~ 
A 
P 



Au 11.64 GeV 



363L10 
1.234±0.126 
133.7L9.9 
23.7L2.9 
3.76L0.47 
20.34L2.74 
0.0185L0.002 



364.1 
1.20 
131.6 
25.77 
3.725 
17.30 
0.0185 



362.9 
1.25 
123.3 
28.55 
3.824 
19.05 
0.0183 



Pb+Pb 404 GeV 

Nb 
A/ A 

A/p (Pb+Au) 

7T _ 
7T+ 

K+ 
K~ 
A 



1121 
I is] 

lis] 

1121 

in] 

2C] 
21] 



349±5 
0.025 ±0.0023 
0.22 ±0.05 
312 ±15 
282 ±15 

56 ±3 
17.8±0.9 
0.71 ±0.07 



351.6 
0.0258 
0.212 
306.2 
276.2 
57.4 
18.2 
0.727 



352.8 
0.0233 
0.232 
264.4 
239.0 
63.2 
19.4 
0.71 



Pb±Pb 804 GeV 

N B 
ir~ 

7T+ 

K+ 
K~ 
A 
A 



349±5 
445 ±22 
414 ±22 
79 ±5 
29 ±2 
47.4 ±3.7 
2.26 ±0.1 



351.6 
403.7 
375.6 
83.0 
33.4 
34.5 
2.26 



352.8 
366.7 
341.1 
87.5 
35.4 
35.7 
2.26 



1 25] 

m 
m 
m 



1 22} 



interaction cross sections are not important here leading 
to a homogeneous strangeness distribution over the whole 
reaction volume. Finally, the SM indicates that the net 
baryon density of the fireball formed in A±A collisions 
at 11.64 GeV and 404 GeV is about 10% below the 
normal nuclear density po = 0.17 fm~ 3 . At these ener- 
gies hydrodynamic model calculations with and without 
the QGP formation |3(| and microscopic model simula- 
tions |2^, [nj favor, however, the formation of a much 
denser matter with central baryon densities of p rs 2-3 po 
at chemical freeze-out, which should take place around 
t = 7-10 fm/c. Thus, the TSM predictions of a small 
central source with baryon densities of ps — 2.5 po quan- 
titatively agree with these estimates. The volume of the 
central source varies from 300 fm 3 at AGS to 650 fm 3 at 
404 GeV H^]. Therefore, to probe this zone one has to 
study hadron abundances and ratios in a quite narrow 
midrapidity window. 

The baryon fraction in the hadrons yield is decreasing 
from 64.5% (59.3%) in the TSM (SM) at AGS to 13.7% 
at RHIC, while the antibaryon fraction rises from zero 
to 9.2%. Figure || depicts the hadronic densities at the 
chemical freeze-out stage in the central zone of heavy- 
ion collisions at ^/s = 1304 GeV and_at 404 GeV. The 
most abundant species are tt, K, p, K, 77, and u> in the 
mesonic sector, and N, A, X, and A in the baryonic 
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TABLE II: Temperature T, baryon chemical potential pb, strangeness chemical potential ps, net baryon density pb, and net 
strangeness density ps obtained from the statistical model fit to experimental data on A+A collisions at 11.6,4 GeV, 40,4 GeV, 
and 80,4 GeV. Of each three numbers the upper one corresponds to the single-source model, the middle number to the central 
core, and the lower one to the halo in the two-source model. 
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123±5 


558±15 


122±4 


0.149±0.01 





Au+Au, 11.6,4 GeV 


141.4±5 


564±15 


142±4 


0.380±0.03 


-0.008±1 x 10" 3 




100.6±7 


558±17 


94±3 


0.004±3 x 10~ 4 


0.0004±5 x 10 -5 




148±5 


367±14 


84±3 


0.142±0.01 





Pb+Pb, 40,4 GeV 


166±5 


413±15 


117±5 


0.430±0.03 


-0.0068±1 x 10~ 3 




103±7 


352±15 


35±3 


0.007±5 x 10" 4 


0.0004±5 x 10" 5 




155±5 


284±15 


66±2 


0.120±0.01 





Pb+Pb, 80,4 GeV 


162. 5±5 


296±15 


74±2 


0.190±0.015 


-0.003±4 x 10~ 4 




98.3±7 


313±15 


28.5±2 


0.003±3 x 10~ 4 


0.0005±5 x 10" 5 



one. Resonances are playing an important role in the 
hadronic spectra at both energies: the average baryon 
and meson masses are 1.22(1.38) GeV and 432(585) MeV, 
respectively, at 40(Vs = 130)^1 GeV. 
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FIG. 3: Densities of (a) baryons and (b) mesons in the central 
zone of A+A collisions at */s = 130 A GeV (stars) and at 40,4 
GeV (boxes) at chemical freeze-out (without feeding). 

The reason for the dropping of the average hadron 
mass and the energy per hadron in Pb+Pb at SPS and 
their subsequent rise in Au+Au collisions at RHIC is as 
follows: If the matter is dominated by baryons, then, de- 
spite the lower temperature at chemical freeze-out, the 
ratio (E)/(N) is larger than 1 GeV. Compared to lower 
energies, the average baryon and meson masses at SPS 
are rising to (M B ) = 1.29 GeV and (M M ) = 470 
MeV, but the baryon+antibaryon fraction drops to 23% 
of the total amount of hadrons. Therefore, the aver- 



age mass of a hadron at chemical freeze-out decreases, 
as was first shown in [|. At RHIC energies the com- 
bined yield of baryons and antibaryons is essentially the 
same as that at SPS, but the average meson and baryon 
masses and kinetic energies increase due to the rise of 
the chemical freeze-out temperature. This leads to the 
rise of the energy per hadron. Note, that if the freeze- 
out temperature increases insignificantly with the rise of 
energy to y/s = 2004 GeV, the combined B + B yield 
remains pretty stable when the baryon chemical poten- 
tial pi b drops to 29 MeV, as predicted in 0, and even 
to zero. Although in the later case the net baryon den- 
sity is essentially zero, the baryon+antibaryon yield still 
amounts 22.8% of the total number of hadrons. Thus, 
the average energy per hadron and the average hadron 
mass should saturate at values (E)/(N) = 1.1 GeV and 
(M) = 770 MeV. 

In summary, average masses (M) and average energies 
per hadron (E) j (N) at chemical freeze-out have been an- 
alyzed for heavy-ion collisions at energies spanning from 
AGS to RHIC. It is found that the rise of the (E)/(N) 
ratio from 1 GeV to 1.1 GeV at RHIC is caused by the for- 
mation of a hot resonance-rich substance, in which the 
fraction of antibaryons becomes significant. Despite of 
the low net baryon density, the combined relative yield 
of baryons and antibaryons is comparable to that at SPS, 
while the average masses of both baryons and mesons are 
increasing due to a higher freeze-out temperature. If the 
temperature of chemical freeze-out will not exceed the 
limit of 180 MeV with rising incident energy, the satura- 
tion values are (E)/(N) = 1.1 GeV and (M) = 770 MeV 
even at zero baryo-chemical potential. A possible viola- 
tion of the scaling (E) /(N) = 1 GeV in A+A collisions at 
energies between 11.64 GeV and 404 GeV has been dis- 
cussed. We show that in baryon-dominated matter with 
baryon densities of 2-3 normal baryon density the energy 
per hadron ratio should rise to about 1.2 GeV. This es- 
timate can be checked by measuring hadron yields and 
ratios in a very narrow range at midrapidity in heavy-ion 
reactions at energies of the planned GSI accelerator. 
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